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1. Introduction 

A satisfactory understanding of long-standing problems in kaon physics, such as the well 
known A/ = 1/2 rule, has so far been elusive. In addition to final-state interactions between the 
two pions, the other possible origins of the A/ = 1/2 rule are "intrinsic" long-distance QCD effects 
at typical energy scales of a few hundred MeV, as well as the decoupling of the charm quark from 
the light quark sector, owing to its large mass of around 1.3 GeV. In refs. ||, |] we outlined a 
strategy to identify a mechanism for the A/ = 1/2 rule, by separately quantifying each of the above 
contributions. Leaving aside final-state interactions, our strategy is implemented by computing 
appropriate hadronic con^elation functions allowing to determine the weak low-energy constants 
(LECs) appearing in the effective chiral theory. Our approach is characterized by the following 
features: 

• The use of overlap fermions [^] in computations of hadronic matrix elements of 4-quark 
operators mediating A5 = 1 transitions. As described in [Q] the mixing with operators of 
lower dimensions usually encountered with Wilson fermions is completely avoided. 

• Matching to ChPT in the so-called £-regime of QCD, where the chiral counting rules imply 
that this step can be performed at NLO without the appearance of unknown LECs. Since 
overlap fermions preserve chiral symmetry the matching can be performed in a conceptually 
easy and clean manner at non-zero lattice spacing. 

• Investigation of the role of the charm quark, keeping it as an active quark in the formulation 
of the effective AS = 1 interaction. This allows to isolate the contributions due to a lai^ge 
mass splitting between nic, ttiu- To this end we start with the (unphysical) situation of a 
mass-degenerate charm quark, = m„ and compute LECs as a function of nic. 

In this note we demonstrate the feasibility of our strategy in the mass-degenerate case, mc = niu = 
md = Tils, where QCD possesses an SU(4)l x SU(4)r chiral symmetry. Since simulations in the 
£-regime aix plagued by large statistical fluctuations ||5|, |6|], we describe in detail how a reliable 
signal can be obtained for 3-point coiTclation functions using "low-mode averaging" (LMA) ||^ ^. 
Furthermore, we discuss the relations between the computed correlation functions and transition 
amplitudes for K nn decays. This requires knowledge of the short-distance renormalization 
factors of 4-quark operators, as well as finite-volume corrections that are computed in ChPT. 

2. A5 = 1 transitions with an active charm quark 

In order to make this note self-contained, we report the basic features of our approach. Ref. [p 
can be consulted for full details. The decay of a neutral kaon into a pair of pions in a state with 
isospin a is described by the transition amplitude 

T{K^^KK\j = iAae'^'', a = 0,2, (2.1) 

where 5a is the scattering phase shift. The experimental observation that the amplitude Aq is sig- 
nificantly larger than A2, i.e. A0/A2 = 22.1, is called the A/ = 1/2 rule. Our task is the computation 
of correlation functions involving local operators, which can be linked to the amplitudes Aq and A2. 
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The relevant local operators ai^e obtained via the operator product expansion of the A5 = 1 
effective weak interaction. For two generations, the effective weak Hamiltonian with an active 
charm quark reads 

•^. = ^{Vu.rV,,d I {^fGf + ^2G?}, (2.2) 
^"^W a=± 

Qi = [{sy^P-u){uy^P^d) ± {sy^P^d){uy^P^u)^ - (m^c), (2.3) 
Q^ = {ml-ml)[md{sP+d) + m,{sP-d)^, P± = i(l±75). (2.4) 

Since Q2 does not contribute to the physical K ^ nn transition we drop it from now on. Note 
that the operators and transform according to irreducible representations of SU(4)l of 
dimensions 84 and 20, respectively. 

The renormalization and mixing patterns of Qf , Qf derived formally in the continuum theory 
are preserved on the lattice, provided that the lattice Dirac operator D satisfies the Ginsparg-Wilson 
relation [||], and therefore an exact chiral symmetry at finite lattice spacing exists If one 
furthermore replaces i/a by ij/ = (1 — ^aD)^f, the resulting local operators in the lattice theory 
have simple transformation properties under the chiral symmetry. Thus, no mixing with lower- 
dimensional operators can occur [^. 

The amplitudes Aq and A2 can be related to low-energy constants in an effective low-energy 
description of A5 = 1 weak decays. To this end we consider the leading order effective chiral 
Lagrangian 

= \F^Tv [{d^U)d^U^] - ^ITr \uM^e!^l^' +MU''t-'^l^''\ , (2.5) 

where U G SU(4) denotes the Goldstone bosons, 6 is the vacuum angle, and M is the quark mass 
matrix. The LECs F and £ denote the pion decay constant and the chiral condensate in the chiral 
limit. The low-energy counterpart of the A5 = 1 effective weak Hamiltonian is obtained at lowest 
order in the chiral expansion as 

=<^'^^ = :|^(K.)*Krf I 8^{muud-msccd}, (2.6) 

where operators containing M have been neglected, and 

W^UprS = i^'^ {Ud^U^)^^ {Ud^U^)gp + projections onto irreps. of dim. 84, 20. (2.7) 

The expression which links the LECS and to the ratio of amplitudes Aq /A2 at leading order 
in ChPT then reads 



A2 V2V2 2g+y 

Finally, the LECs can be determined by matching suitable coiTclation functions in ChPT and 
QCD. This leads to 

g^ ^_ k-iMw/A) Z-(go) 

gi k+iMw/A) Z+{go) Cf 

Here, the chiral correction factor H is obtained as a ratio of correlation functions of [^^] computed 
in ChPT, and Cj*" are specified in eq. (3.2) below. On the RHS the short-distance corrections include 
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the Wilson coefficients and the renormalization factors Z*, which relate the unrenormalized 
operator (2i^)baie> considered at bare coupling gQ, to the renormalization group invariant operator 
via 

(Gf)RGI = Z*Uo)(ef)bare. (2.10) 

In the following sections we describe the evaluation of the correlation functions, the chiral connec- 
tion and renormalization factors. 



3. Lattice set-up in the SU(4)-symmetric case 

Since preserving chiral symmetry is an essential feature of our setup, the computation of the 



correlation functions in Eqs. ( [3.1|J3.2| ) is performed in an overlap lattice regularization. In order 
to match QCD to its effective low-energy description in the SU(4)-symmetric case, we start by 
defining suitable two- and three-point correlation functions of left-handed currents and four-quark 
operators in QCD, namely' 

C(xo) =£([/oW]a/3[-/o(0)]/3„> , (3.1) 

X 

Ct{xo,yo) = Y.{[Ux)]du (0) [My)]us) , (3.2) 

where the non-singlet left-handed current is defined through 

[Jfii^)]ap = {WayiiP-Wp){x) , (3.3) 
05, j3 are generic flavour indices, and the replacement y ~^ ¥ has been performed in the four- 



quark operators of Eq. ([2.3|). Recall that in the SU(4)-symmetric limit the three-point functions Cj 
receive contributions from "figure-8" diagrams only, since "eye" diagrams exactly cancel due to 
the antisymmetrization under {u ^ c). It is also useful to define the following ratios of correlation 
functions, which will enter the determination of low-energy constants: 

p ^ \ Ct{xo,yo) „ f . Cf(xo,jo) 

C{xo)C{yo) C{xo)C{yo) 
C+{xo,yQ) 

R84/2o{xo,yo) = -p;—, r • (3.5) 

{xQ,yo) 

Far enough from the location of the source operators, all these ratios are expected to exhibit 
plateaux that can be fitted to a constant, which can then be used as input in the matching procedure 
to Chiral Perturbation Theory. 

At low quark masses the numerical computation of correlation functions is usually hampered 
by the presence of large statistical fluctuations. The latter can be understood by considering the 
expression of the quark propagator in terms of eigenmodes of the Neuberger-Dirac operator D, viz 

S{x,y) = -l^ , (3.6) 



The use of left-handed currents, as explained in ||10l Ip, is particularly convenient for technical reasons 
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with At = (1 — \am)Xk and Dr\k = ^k^k■ In the regime m < (ry)^^ which allows a matching to 
Chiral Perturbation Theory in the £-regime, the low-lying spectrum of D,,, = (1 — ^dm)D-\-m is 
discrete with AA 1 /£V, and sizeable contributions to coiTclation functions come from a few low 
modes. Large statistical fluctuations can be traced back to "bumpy" structures in the wavefunctions 
of these modes [12, ^. 



In order to treat this problem we use low-mode averaging (LMA) introduced in The 
technique proceeds by treating expUcitly the contribution to left-handed quark propagators coming 
from a few lowest-lying modes of D. To be specific, we split propagators as 

S{x,y) = > +S {x,y) , (3.7) 

where S^' is the propagator in the orthogonal complement of the subspace spanned by the ?iiow 
lowest modes, et = PaUt + P-aDPaUk^ — <7 being the chirality where D possesses zero modes (if 
any), and is an approximate eigenmode of D\fim'. 

PaDlDjnPaUk = UkUk + rt , {uk,ri)=0 Mk,l. (3.8) 



After inserting the RHS of Eq. (3.7) in the expressions for the correlation functions C and Cj*^, they 
can be split as 

C = d^ + C''^ +d''\ (3.9) 

^± ^±\UU _|_^±;////i _^^±\Uhh j^^±\lhhh _^^±;hhhh ^ ^q-j 

where / and h denote the number of "light" and "heavy" parts of the quark propagator, respectively. 
Since the "Ught" part of S{x,y) is available by construction 'ix,y, it is possible to exploit transla- 
tional invariance to sample the all-/ contributions over many different source points. Furthermore, 
as explained in [^, by performing n\o^ additional inversions of the Dirac operator it is also possible 
to extend this to the mixed contribution C'''. It is easy to check that the same applies to the hill con- 
tribution to Cj*^, as well as to part of the hhll one. As already shown by the exploratory study in [pl|], 
the application of this technique with ?iiow ~ 20 suffices to obtain a signal for three-point functions 
at values of the quark mass of interest in view of matching to £-regime Chiral Perturbation Theory 
results. 

4. Correlation functions in the e and p-regimes 

Our simulation parameters are summarized in Table |l|. The simulations for lattice A are those 
reported in [[I]], while lattices B and C ai^e new results. The statistics of lattice C is currently being 
increased. The results quoted for lattices B and C have to be considered preliminary. 

Our main aim is to fit to constants the plateaux in the ratios of coiTclation functions in Eqs. 



( p.4|j3.5| ). For quark masses in the /7-regime the procedure is straightforward, and our statistics 
allows quite precise results for the different ratios. An example for /?84/20 in th^ p-regime is shown 
in Fig. |l]. It also shows the effect of LMA on correlation functions at typical p-regime masses. 

In the £-regime topology plays a special role [pj|], and coiTclation functions are given within 
fixed topological sectors. Therefore we proceed by computing the quantities of interest at fixed 
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Lattice 


/3 


L/a 


T/a 


'^low 


L[fm] 


am 


# cfgs 


A 


5.8485 


12 


30 


5 


1.49 


0.040,0.053,0.066,0.078,0.092 


638 


B 


5.8485 


12 


32 


20 


1.49 


ft OOJ, ft 005, 0.007,0.040 


681 


C 


5.8485 


16 


32 


20 


1.99 


ft 002, ft 00i,0.020,0.030,0.040,0.060 


~350 



Table 1: Simulation parameters for the runs discussed in the text. The mass values in italics are those 
corresponding to the e-regime. The statistics indicated for lattice C refers to the number of configurations 
for masses in the e-regime; for p-regime masses the statistics is roughly half of the indicated figure. 




0.2 0.4 0.6 0.8 1 



0.0 0.2 0.4 

^54/20 



Figure 1: Left: The ratio 7?g4 for am = 0.020, lattice C. The points indicated by circles (triangles) have been 
obtained with(out) LMA. Right: Weighted average over | v| (solid band) of the ratio /?84/20 for = 0.003, 
lattice C. 



value |v| of the (absolute value of the) topological charge, and then perform a weighted average 
over |v|. In order to have large enough statistics within each sector, and taking into account the 
expected distribution of topological chai^ges, we impose a bound on the largest value of | v| entering 
the average (| v| < 8 on lattice B and |v| < 10 on lattice C). Furthermore, following the observation 
that the signal-to-noise ratio in the sectors with lowest |v| is poor,^ for the largest volume (lattice 
C) we also impose a lower bound |v| > 2. This procedure is illustrated in Fig. |TJ In the £-regime 
we find no signal at all for the relevant observables if LMA is not implemented. Indeed, to our 
knowledge, these are the first results for three-point functions obtained at quark masses in the 
e-regime. 

Our most interesting results ai^e those for the mass dependence of the different ratios. They 
are summarized in Fig. ^, where we put together the /^-regime results for our three lattices and 
the £-regime results in our larger volume. Three features worth mentioning are: Firstly, the mass 
dependence is remai^kably smooth. In particular, there is no strong mass dependence of /?84 for 
very small quai^k masses. Notice, however, that finite volume connections have to be taken into 
account for the £-regime points (see below). Secondly, our results point towards a moderate volume 
dependence at quark masses corresponding to pseudoscalar meson masses around or below the 

^This can be interpreted as a consequence of the presence of very small eigenvalues of D, which in turn induce large 
statistical fluctuations even after LMA. 



344/6 



On the determination of low-energy constants for AS =1 transitions 



Carlos Pena and Hartmut Wittig 




0.00 



1.0 

0.8 

0.6 

0.4 - 

0.2 F* 

0.0 
0.00 



0.08 




0.0 



0.00 



0.02 



0.04 0.06 
am 



0.08 



0.02 



0.04 
am 



0.06 



0.08 



Figure 2: Mass dependence of the ratios Rgn, R20 and Rs4/2q for lattices A (crosses), B (empty circles) and 
C (full circles). Results for lattice C are preliminary. The points at am = 0.040 have been slightly displaced 
to improve visibility. 



kaon mass. As far as Rg4 is concerned, this was akeady observed in [|T^]. Thirdly, the direct 
comparison of the different results at am = 0.040 shows that the effect of LMA with an adequate 
number of low modes is far from negligible even at moderately large values of the pseudoscalar 
meson mass. 



5. Renormalization factors for 4-quark operators 



The renormalization factors Z*(go) of eq. ( 2.10 ) are scale and scheme independent. For a 
particular renormalization scheme X they can be decomposed according to 



) = c| ( At / A) Z| (^0 , a At ) , 



(5.1) 



where jLt denotes the renormalization scale, and the coefficients are given by 



r±fl/{2ho) 



exp 



dg 



T^(g) y±;o 

Pig) bog 



(5.2) 



The anomalous dimensions are known in perturbation theory to two loops for several schemes. 
For discretizations based on the Neuberger-Dirac operator, the renormalization factors Zx{go,aiJ.) 
have been computed for X = RI/MOM in perturbation theory at one loop in ref. [Q]. Thus, the 
perturbative renormalization of suitable ratios of 4-quark operators defined for overlap fermions 



344/7 



On the determination of low-energy constants for A5 = 1 transitions 



Carlos Pena and Hartmut Wittig 



and the RI/MOM scheme is 



4;Rlteo) 



l + Y^{121n(4Ata) 



41n(4/ifl') 



2iBs-By)} + 0igt,) 



1671^ t 



2i;Rlteo) 



1 + 



16712 



81n(4/xa) 



3(5s- 



■fiv)!>+ofeo) 



(5.3) 



(5.4) 



(5.5) 



The coefficients 8$ and By are listed in Table 1 of [Q]. In addition to 4-quark operators, we have also 
considered the renormalization of the axial current. Eqs. (5.4) and (^jsb 



then serve to renormalize 

the corresponding B-parameters of the operators Qf. The RGI matrix elements aie obtained by 
combining the above renormalization factors with the coefficients c^. 

Perturbation theory in the bare coupling gg is known to have bad convergence properties. 



The aim of "mean-field improvement" [15| is to factor out unphysical tadpole contributions in the 
perturbative expansion, by a rescaling of the link variable, U^{x) U^{x)/uo. For the Neuberger- 
Dirac operator defined by 



^(1 



l+s, \s\ < 1, 



1 + 5 — aDy, 



(5.6) 



the corresponding rescaling of the quark field is given by xjf ^ \/{p/p)¥: P = {P — 4)mo + 4. 
For the renormalization factor = 1 +§0^^' + 0(go) of an n-quark operator ff„, the mean-field 
improved version reads 



n/2 



n p 
2~ 



,(1) 



(5.7) 



where uq 



l + «i>) 



When applied to our set of operators, it is immediately clear that the 

(1) 



contributions from the prefactor (p/p) as well as those proportional to Uq drop out in ratios like 



/Z^ and Z'^/Z^. Mean-field improvement of the expressions in eqs. (53)-(5^) is thus simply 
accomplished by replacing the bare coupling by the "continuum-Uke" coupling g-^. 

For a reliable determination of operator matrix elements, the use of non-perturbative estimates 
for renormalization factors is to be preferred. The Schrodinger functional (SF) offers a general 
framework for non-perturbative renormalization of QCD at all scales [p^]. However, the con- 
struction of SF boundary conditions consistent with the Ginsparg-Wilson relation is quite involved 



[|17p. In ref. 018| ] it was therefore proposed to introduce an intermediate Wilson-type regulariza- 
tion which drops out in the final result. As an example we now discuss the renormalization factor 
Z'^ lz\, which is required for the B-parameter Bk- The desired factor relates the B-pai^ameter 
S^(go) computed using overlap fermions, to its RGI counterpart Bk- After introducing an inter- 
mediate regularization based, for instance, on twisted mass QCD [[l9|], it can be written as 



Bk 



Bk Bf{g',) 



Z^ 

'^A;ov 



B7i 



Bfig'o) B°^{go 



lim 



■Bf{ 



B7i 



(5.8) 



where the superscripts "ov" and "tm" on the unrenormalized B-parameters refer to overlap and 
twisted mass fermions, respectively. The key observation is that the expression in square brackets 
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is nothing but Bk in the continuum limit, which, for instance, has been computed by the ALPHA 
Collaboration in quenched QCD [14]. Denoting the result by 5^^™^, the renormalization factor 



in eq. (5^) is B^^^^ / B'^^ [go) . Of course, in this way one cannot predict Bk any more, since its 
value is used to formulate the renormalization condition. However, the procedure can be used to 
determine the value of Bk in the chiral limit, B^, in units of B^^^^: 

Bi = X 4- Uo) + 0{a^) = ^r"^ X + o{a\ (5.9) 

2a;ov \So) 

Note that Bf'°^(go) can be obtained from a suitable ratio of correlators computed in the e-regime, 
in conjunction with the appropriate chiral correction factor. 

We now discuss some numerical examples 
for perturbative and non-perturbative estimates 
of renormalization factors. In our simulations 
we use j3 = 6/gl = 5.8485. For = 2GeV and 
A = 23 8 Me V [^0|], the perturbative expressions 

for the coefficients c± yield CRi(At/A) =0.6259 

J +t /\ \ 1 noc ivT 4^ I. 4^- Table 2: Perturbative and non-perturbative estimates 

and u/A = 1.2735. Non-perturbative es- ^ 

, . for RGI renormalization factors at j3 =5.8485. 
timates for the B-parameters computed at the 

physical kaon mass in the continuum limit of quenched QCD are provided by the ALPHA col- 
laboration [14]. The results for ratios of renormalization factors are listed in Table ^ 







bare PT. 


MFI PT. 


non-pert. 


z 


7Z+ 


0.525 


0.582 


0.58(8) 


Z^ 


74 


1.242 


1.193 


1.20(8) 


Z 


74 


0.657 


0.705 


0.73(8) 



The entries in the table show that non-perturbative estimates for renormalization factors are 
remarkably close to perturbative ones. Indeed, even the differences between perturbative estimates 
evaluated in "bai^e" or "mean-field improved" perturbation theory are small, presumably since ra- 
tios of operators are considered here. This is in stark contrast to the situation encountered for sim- 
ple quark bilinears, for which the deviations between perturbative and non-perturbative estimates 
amount to about 30% at similar- values of the bare coupling [^. 



6. Chiral corrections 

Our strategy of determining the LECs of the AS' = 1 weak interactions requires that the kine- 
matical range where ChPT is applicable must be accessible to lattice simulations of QCD. The so- 
called e-expansion [ p2| ] represents a systematic low-energy description of QCD in a finite volume 
for arbitrarily small quark masses. It is characterized kinematically by the conditions mLV 0(1), 
Fy^L » 1, where V = is the four- volume, and m is the quark mass. These conditions lead to 
different chiral counting rules compared with the more commonly known /^-regime. In particular, 
since the inverse box size counts as one unit of momentum L^^ 0{e), one infers m ~ 0(£^) and 
hence mj^ ~ 0(£^). For the effective Hamiltonian J^^^^ of eq. (2.6) this in turn implies that no 



additional interaction terms are generated at 0(e ). In other words, the £-regime allows for a NLO 



matching of lattice data to ChPT without the appearance of additional, unknown LECs [23]. 



We can now work out the chiral correction factor H in eq. (2.9). To this end we define corre- 
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lation functions of the left-handed axial current in complete analogy with the fundamental theory: 

^"\xo) = J d3^(^o"WX(0)) > A" = -2F\n,p (6-1) 

l'^^i^o,yo)Cy5 = I / ^'y{/S{^)Wxma^ys/^{y)) , (6.2) 



Choosing a diagonal quark mass matrix and flavour matrices T",T'^ as in eq. (D.6) of one 
defines the chiral correction factor H by 

■^j {xo,yo) 

H = ^— = l-2R[xo,yQ). (6.3) 

For later use we also consider the chiral corrections for B-parameters, i.e. 

K^^2 =l + aR{xo,yo), o = ±. (6.4) 

^(xo)^(3'o) 

Explicit expressions are listed in section 5.3 of ref. ||l|]. In Fig. 2 of [[I]] the quantity R is plotted as a 
function of the box size for several lattice geometries. It clearly demonstrates that chiral corrections 
are reasonably small for box sizes L > 1.5fm and lattice geometries with T /L <2. 

7. Synthesis, conclusions and outlook 

We can now combine our results for ratios of coiTclation functions with the appropriate renor- 
malization and NLO chiral correction factors. We expect that the latter are best controlled for our 
dataset "C", for which T /L = 2 and L = 2fm. The link between the ratio gi /g^ and the con^elation 
functions is given in eq. (^^), while individual values for g'f are obtained from 



zHgo) 



giK+ = k+{Mw/A)-^ ■ (^o), (7.1) 

and similarly for g^ . The LECs gf are then related to the amplitudes Aq, A2 via LO ChPT. 

Our preliminary results for these amplitudes in the SU(4)-symmetric theory indicate a severe 
mismatch with experiment: roughly speaking, our value for Aq is too small by a factor 2, while A2 
comes out a factor 2 too large. This produces an estimate for A0/A2 which is four times smaller 
than the one expected from the experimentally observed A/ = 1/2 rule. On the other hand, this is 
a factor 4 larger than the naive large-A'^c limit, and does thus move in the right direction compared 
with this case. 

However, it would be premature to conclude that the A/ = 1/2 rule is generated by the decou- 
phng of the charm quark, since the amplitude A2 is insensitive to the charm mass, yet its experimen- 
tal value is not reproduced either in our calculation. Other possibilities for the observed mismatch 
are uncontrolled finite- volume corrections, quenching effects, or even the breakdown of LO ChPT 
when relating the LECs to the transition amplitudes. Our future work will thus concentrate on 
corroborating our results in the e-regime, as well as incorporating the effects of a non-degenerate 
charm quark mass. In this context we shall investigate alternative choices of correlators, which are 



saturated with zero modes [24|. 



Our calculations were performed on PC clusters at DESY Hamburg, CILEA and the Univer- 
sity of Valencia, as well as on the IBM Regatta at FZ Jiilich. We thank all these institutions and 
the University of Milano-Bicocca for their support. 
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